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We present a specific class of models for an infrared-finite analytic QCD coupling, such that at 
large space-like energy scales the coupling differs from the perturbative one by less than any inverse 
power of the energy scale. This condition is motivated by the ITEP Operator Product Expansion 
philosophy. Allowed by the ambiguity in the analytization of the perturbative coupling, the proposed 
class of couplings has three parameters. In the intermediate energy region, the proposed coupling 
has low loop-level and renormalization scheme dependence. The present modification of perturbative 
QCD must be considered as a phenomenological attempt, with the aim of enlarging the applicability 
range of the theory of the strong interactions at low energies. 
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I. INTRODUCTION 



The construction of an analytic coupling Ai{Q^), a coupling analytic in the Q^-plane excluding the negative 
(time-like) semiaxis, is a manner to provide the evaluated observables with the analytic properties required by 
causality. Perturbative QCD (pQCD), based on a truncated power series of the /3- function, yields a coupUng 
flpt(Q^) = cits{Q'^)/n which contains (Landau) singularities in the low energy space-like regime (0 < < A^). These 
singularities are problematic as they result in nonphysical singularities of the space- like QCD observables 2?(Q^) when 
the latter are evaluated as truncated perturbation series (TPS) in the pQCD coupling apt(/.t^) (with renormalization 
scale = kQ'^, k ~ 1). The most direct construction of an analytic coupling was performed by Shirkov and 
Solovtsov [l| (minimal analytic coupling - MA; for applications, see 0,13, 01): who kept the discontinuities of apt(Q^) 
unchanged on the time-like axis, but removed them from the space- like axis.^ Other models for an analytic coupling 
have been proposed afterwards [a, [l^, [HI, > which change in general the low energy behavior with respect to the 
MA coupling. For a review of various models, see [Tsj . 

At high energies > A?, all these coupHngs differ from apt(Q^) by '-^ (A^/Q^). In such models the power 
suppressed terms in space-like QCD observables would come, at least partly, from the ultraviolet (UV) regime. This 
would contradict the philosophy of the ITEP Operator Product Expansion (ITEP-OPE) approach [IJ]. If we want 
the ITEP-OPE approach to survive in analytic QCD models, then 5Ai{Q'^) = -4i(Q^) - apt((5^) at large (> A^) 
must fall faster than any inverse power of . We will show that, within the context of one-chain resummations, even 
in the case when 5Ai{Q'^) at ^ A^ is a very suppressed power correction 5Ai{Q'^) ^ (A^/Q^)'^"''^ (A:max ^ 1), 
the UV regime still contributes to the space-like observable 'D{Q'^) an appreciable power ~ (A^/Q^)". Here z — n 
is the location of the leading IR renormalon of V in the Borel plane. On the other hand, such term has the same 
power-behavior as the leading OPE term of the observable, which is of IR nature. This case would thus contradict 
the ITEP-OPE philosophy. In the present paper we shall study this case in detail and we shall see that the definition 
of UV and IR contributions in the renormalon-like resummation is a crucial point. 

In Ref. 15], a coupling with \5Ai{Q'^)\ < ©((A^/Q^)"^) is constructed, by adding two specific power-suppressed terms 
to the MA-coupling. In Ref. [ly|, a class of analytic couplings is obtained from modified /3-functions. The obtained 
SAi{Q^) is at large smaller than any inverse power of Q^. Both types of couplings are divergent at ^ 0. 

In this work we construct a class of infrared-finite analytic couplings Ai{Q^). The proposed couplings deviate 
from the perturbative one by terms exponentially suppressed at high Q^, thus they differ from UptiQ^) by less 
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than any inverse power of Q^. While in the UV region Ai{Q'^) mimics the perturbative coupHng, the IR behavior 
of Ai{Q^) is modeled by means of a set of parameters (three). In the intermediate energy region, the proposed 
coupling has low loop-level and RS dependence. The reason why the coupling is modeled lies, on one hand, 
in the fact that there is not a unique way to analytize the perturbative coupling and, on the other hand, in 
the possibility of enlarging the applicability range of pQCD by maximizing the description power of the theory at 
low energies. This must be considered a phenomenological attempt where the relevant energy region is Q ~ 1 — 2 GeV. 

In Section 2, we present the long-distance and the short-distance sources of the power corrections (A^/Q^)" at the 
level of the one-chain resummations (the leading skeleton, or equivalently, the large-/9o approximation). In Section 3, 
we present the construction of a class of such infrared-finite space-like analytic couplings Ai{Q^) which differ from 
apt(Q^) by less than any power (A^/Q^)" for large Q^. We also present the associated Minkowskian coupling. In 
Section 4, we discuss the loop-level dependence and renormalization scheme (RS) fixing of low-energy observables. 
We also review the method of evaluation of observables developed in our previous works [111, [l2| for general analytic 
QCD models, and see how it can be applied in the case of the here constructed analytic coupling. Section 5 contains 
the conclusions of this work. 



In this Section we consider different contributions to the one-chain resummation term I?rcs(Q^)7 paying special 
attention to their IR or UV character. We are mainly interested in power-behaved contributions. 

The one-chain resummation term is given by 



where F^{t) is the observable-dependent Euclidean characteristic function and Ai{Q'^) the observable- independent 
running coupling, normalized as a(Q^)/7r. (The constant C is the one-loop vacuum polarization renormalization 
constant. It has to be included in order to make I'rcs(Q^) C-independent. In the following we use the V-scheme, i.e. 
C = 0, unless otherwise noted.) We consider a running coupling, Ai{Q'^), analytic in the whole complex plane 
excluding the Minkowskian semiaxis, i.e. a so called analytic coupling. The analytic coupling must merge with the 
perturbative one in the UV limit. In the UV expansion of Ai{Q'^) one finds, in general, in addition to the usual 
logarithmically varying terms of the perturbative treatment, also power-behaved terms of the form (l/Q^)", with 
n — 1,2, . . .. We call them running coupling power terms. The ^-integration in Eq. ^ is an Euclidean momentum 
integration, where the running coupling is evaluated at the flowing momentum scale, t^/^Q. This momentum is 
identified with the momentum passing through a photon/gluon propagator. In QED there is a correspondence 
between vacuum polarization and charge renormalization which allows the identification. In QCD there is no such a 
correspondence and the identification is naive. 

We aim to identify IR and UV contributions to I?rcs(Q^) by comparing the momentum flowing through the boson 
propagator, t^^^Q, with the scale where QCD becomes non-perturbative. For simplicity, we make a sharp division 
between IR and UV regions. The contributions coming from the integration region where tQ^ < kA^ are considered 
to be IR, while those coming from tQ^ > kA? are considered to be UV. The constant k is chosen to be bigger than 1, 
and of the order or a couple of orders of magnitude bigger than \? Note that these IR and UV regime identifications 
are different from those where the IR and UV regimes are identified as i < 1 and f > 1, respectively. If ^ A^, the 
IR and UV regions are similar in both IR/UV identification criteria, but differ for » . 

In order to be concrete, we take as an example the simplest analytic coupling: the one-loop minimal analytic 
couphng [l|: 



^ For the treatment below, the specific value of ft is irrelevant, it can be chosen to be 2, 10 or 100; of importance is its Q^-independence. 
The use of a sharp division between IR and UV regions is also not decisive. The point is that, using this criterium, the minimal value 
of t in the one-chain resummation integral Eq. JTJ from which the contribution is considered to be UV is Q^-dependent. 



II. POWER CORRECTIONS FROM ONE-CHAIN RESUMMATIONS 




(1) 
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4""^) (g2) _ ( ^^^^^ Q2 _ aO ' 

where the first term is the perturbative one-loop running coupUng, apt{Q^) = apt(Q^)/7r = (/3o logQ^/A^)"-'^, and 
the second term is introduced in order to make Af^''^\Q^) analytic in the mentioned region without changing the 
discontinuity along the Minkowskian semiaxis. We inspect now the contributions to 2?res(Q^) from both terms 
separately. 

A. Perturbative Coupling Contributions 

The perturbative coupling contribution to I?ros(Q^) is 

We shall consider different aspects of the latter expression. First we obtain its contributions to the usual perturbation 
series, expanding the log^^ inside the integral. Then, we give the usual definition of the integral ambiguity, a 
reflection of the fact that the integral in Eq. ^ is not a well defined quantity. Finally, we consider the contribution 
from the IR integration region. 

In order to obtain the part of the perturbation series which is contained in Eq. ([3]), we expand the running coupling 
o-ptitQ^) around a chosen renormalization scale Q^: 

oo 

apt(iQ') - 5^(-/3ologiA*)"<t+'(i*Q')- (4) 
Replacing in Eq. ([3]) and exchanging the order of integration and summation, we obtain 



where 



UU)^ / -F|(t)(-logtA,)"- (6) 

The series obtained in Eq. ^ is the large- /3o expansion of which consists of all the terms of the perturbation 

series of the form /3Japt+^ If we had considered in Eq. ^ another perturbative running coupling, solution of a 
higher order renormalization group equation, we would have obtained some further terms (beyond the large-/3o) of 
the perturbation series of the observable 'D{Cf'). 

The integral in Eq. ^ is not well defined due to the simple pole ai t — jQ"^. Clearly, replacing apt((5^) by 
an analytic coupling AiiCf'") one gets a well defined quantity as our starting point, Eq. ([1]), is. In order to give 
a meaning to I?Jgg((3^) one chooses a modification of the integration contour, relevant only at the pole. The most 
used prescription is to take the principal value of Eq. (O; the procedure, however, yields a non-analytic but piecewise 



analytic contribution, where the non-analytic term is purely imaginary [17| . Related to the prescription choice there is 
an ambiguity, the standard definition of the ambiguity is: the difference between Eq. ([3]) with the integration contour 
slightly above the real axis and Eq. ^ with the integration contour slightly below the real axis. As a result we get: 

AmW'\v\,,{Q')\ = F%{k^lQ\ (7) 
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In the following we consider the small t behavior of (t) to be given by 

F^{t)^kt", fort<l. (8) 

Thus, for > one gets 

= + (9) 

The ambiguity Eq. ([7]) contains the contributions of all large- /?o IR renormalons. These poles lie on the Borel plane 
at u = n where n is a positive integer [Tsj . For the relation between the representations of the large- /3o contribution, 
using a one-chain resummation or using a Borel integral, see [l^. In the present formulation, the ambiguity ([7]) 
depends on the characteristic function and on the (non)analyticity properties of the perturbative running coupling 
- Landau singularity (pole at tQ'^ = A) [^oj. For inclusive observables, such as sum rules, the ambiguity of the type 
(171) must get canceled by the ambiguity of the corresponding power term in the Operator Product Expansion (OPE) 

The contribution from the IR integration region, t < kA^/Q^, is given by 



For ^ A^ its leading term is ~ (A^/Q^)". The latter integration region includes the position of the running 
coupling pole and Eq. pH)) has the same leading Q^-dependence as the ambiguity. 



B. Contributions from Power Terms 



The difference between the I'rcs(Q^) and the perturbative coupling contribution is given, in the case of the one-loop 
minimal analytical coupling, by 

^ ^.,AQ') - vijQ^) - [ . (11) 

The integrand has a simple pole at the same value as the integrand in Eq. ([3]). The corresponding ambiguity, 
+{2m/ Pq)Fj;,{A? /Q'^), has equal absolute value and opposite sign as the one of I?Jgg(Q^), therefore cancels in 
Pros = 2^Jes + ^rcs- This implies that, if OPE is applied in the analytic approach, the OPE terms must also have 
their ambiguity lifted. One can also see that the IR part of Eq. (fTTj) has the same leading Q^-dependence (A^/Q^)" 
as the IR part of I?J^,(g2), for > A^. 

Now we turn to the aspect we are mainly interested in, the UV contribution: 

dt f . . 1 /■ -K 



Since k > 1, inside the integral we are allowed to expand: 



We consider the contribution to ^^^{Q^) from each term of the latter sum individually 



if ?^^«(i)'(^ 



1\ V A2\ 

^ (14) 
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Then, we ask ourselves what is the Q^-dependence of the latter expression for ^ A^. In order to use the 
approximated expression for F^{t), for t <C 1, we divide the integration interval in two parts: 

/ = / + / , (15) 

choosing (fixing) t < 1 such that in the first integral the approximation for F^{t) given in Eq. ([5]) is valid. We see 
that the second integration interval generates a power correction term of the form (A^/Q^)*, while the contribution 
from the first integration interval can be written as: 

1 r dt /1^VA^^ r (AW)Mog(Q W + . . . if z = n, 

Thus, the term ([14]) has a behavior given by 

I r dt /1^VA^^ r(AW)Mog(QVA^) + ... ifz^n. 

We conclude that terms of the type (l/Q^)' from the high momentum expansion of an analytic coupling Ai{Q'^)^ 
introduce in the UV integration region of the one-chain resummation integral, Eq. ([1]), observable power corrections 
whose leading term is proportional to {A'^/Q^Y log(QVA^), if i = n, and to (AVQ^)™"^*'"\ if i ^ n, for > A^. 

Returning to expression (fT2|) . i.e. the UV part of Pjjg, the leading power-behaved contribution comes from the 
term with i = 1; and is proportional to A'^/Q^ if n > 2, or to (A-^/Q^) log(Q^/A^) if n = 1. 

Power corrections coming from UV degrees of freedom are in conflict [l^, [23| with the OPE philosophy. Therefore, 
an analytic coupling having a high momentum expansion with corrections to the perturbative running coupling 
falling faster than any negative power of is wished. We address this question in the next Section. 

Before doing this, let us consider a different definition of the ambiguity in the one-chain resummation term I'J(,g(Q^), 
motivated by the fact that the perturbative running coupling apt(Q^) has no physical meaning in the IR region. We 
can define the ambiguity of 2?];gg(Q^) as the contribution of the IR region, i.e. t < kA^/Q^, with apt{Q^) inside the 
integral replaced by a constant qq: 

Amh[Vl,{Q^)]=aJ -PUt). (18) 

Jo ^ 

First, we note that the ambiguity is now a real quantity, instead of a pure imaginary one as it is conventionally 
defined. Furthermore, for 3> A^, the usual ambiguity and Eq. (fT8|) have the same leading Q^-dependence 
~ (A^/Q^)". The ambiguity defined as in Eq. (fTS]) is independent of the properties of a^tiCp), in particular, is inde- 
pendent of the loop order the perturbative running coupling is solution of, and independent of its analyticity properties. 

In order to have an idea of the order of magnitude of the ambiguity Eq. we take as an example the Adler 

function, Q = 1.8 GeV , Ac=o = 0.9 GeU, and k = 1.5. We get Amb[Vl^^] = 0.8 oq. The ambiguity defined in this 
manner serves as a measure of the effect of modeling the running coupling on one-chain resummation contributions. 

III. EXPONENTIALLY MODIFIED COUPLING 



In this Section a new analytic coupling is proposed. In the 3> A^ region, the coupling converges to the 
perturbative running coupling with deviations from it falling faster than any power of the momentum. The behavior 
of the coupling in the IR is described through few parameters. 
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A. Logarithm Replacement 

We start considering the one loop case. In Section |lT] we saw a possible replacement of the logarithm which yields 
an analytic coupling: 

111 

(19) 



log X log X 1 — a; ' 

where x = / . The latter analytization of the l/(loga;) function is unique if one requires that: (a) it keep the same 
discontinuity along the a; < semiaxis and (b) it vanish for \x\ — > oo (with x complex and not lying on the negative 
semiaxis; and having the dispersion relation shown below). As we discussed above, Eq. (jl9p introduces unwished 
power terms in the large- a; expansion of I?ros(Q^)- In order to modify the latter replacement, we must abandon at 
least one of the previous conditions. It seems unreasonable to give up condition (b). If the coupling does not vanish 
for |a;| — > oo in some region of the complex plane, sum rules and the expression for the Adler function in terms of 
the R ratio would not hold. It is also against the intuition of asymptotic freedom. On the other hand, there is no 
reason to maintain the same discontinuity on the time-like axis as in perturbation theory, hence we relax this condition. 

The following replacement is an example of modification of Eq. (jl9p . that does not introduce power-behaved terms: 

1 1 e^^-^"^") 



log X log X 1 — X 



(20) 



with V > and < a < 1/2. The parameter a must not exceed the value 1/2 because otherwise condition (b) would 
be violated. We define the cut of the a;" function along the negative semiaxis as for the logarithm function. Thus, 
by means of the exponential function, we get a one-loop analytic coupling with corrections to the perturbative one 
falling, in the UV region, faster than any power of x. 

B. TV- loop Coupling 

Based on the previously introduced replacement, Eq. (I^Dl) . we construct an analytic coupling from the iV-loop 
running coupling. We write the latter in the form: 

4t ^ = E E (21) 

n— 1 m— 

where L = logQ^/A^ and knm are functions of the /9-fimction coefficients. The new function A^^\Q'^) should have 
the following properties: 

1. be analytic in the complex plane excluding the negative semiaxis, 

2. differ from a^^\ for 3> A^, by terms falling faster than any power of A^/Q^, 

3. vanish for |a:| oo as a:^^ (with x = Q^/A^ complex and not lying on the negative semiaxis) and have a 
dispersion relation: 

r^pM, (22) 



with p((t) = lm{A^^\-a - ieA"^)), and 
4. suppressed iV-dependence (loop-level dependence) in the IR region. 
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FIG. 1: The functions Eq. l|27ll, and a^^ \ Eq. are plotted as a function of Q, in the MS-scheme, with A = 0.4 GeV, for 

the cases A'^ = 2, 3, and 4. See comments in the text. 

The loop-level and the renormalization scheme (choice of P2, /^a, • ■ •) are perturbative concepts. Demanding the last 
property is necessary for a coupling largely independent of these concepts in the IR region. 

We propose an analytic coupling having the form 

N n~l , „i r 

A'^HQ') = Y.Y. + e-^v^/(x). (23) 

n=l m=0 

The second term is only relevant in the IR region and the first term (double sum) plays, in the UV region, the role 
of the perturbative coupling. Lq and Li are chosen aiming at a low TV-dependence in the IR region. That is achieved 
by suppressing logLi and I/Lq in this region. Generalizing Eq. (|20p we define 

1 1 e^ii^-V^) 

— = - + — g^{x), I., >0, z = 0,l. (24) 

Li L 1 — X 

The functions gi{x) and the constants Vi are chosen as follows. For Lq we require: 

(i) 1/Lq{x) analytic in the complex plane excluding the negative semiaxis. It implies 5o(l) = 1; 

(ii) 1/Lq{x) suppressed (^ 1) near x ~ 1, 

(iii) 1/Lq{x) suppressed near x — 0, and 

(i'^) 9o{x)/x < constant, for \x\ — > 00, in order to fulfill the previous condition 3. 
Similarly, for Li we require: 
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FIG. 2: Same as in Fig.[T] for a different range of energies. 



(i) logLi(a;) analytic in the complex plane excluding the negative semiaxis. It implies gi(l) = 1, 

(ii) logLi(a;) suppressed near x — I, and 

(iii) logLi(2:) suppressed near x = 0. 

We present an analytical coupling fulfilling all previous requirements. We choose: 

2x 

9q{x) = T— — r— 7^ 7, < I/O < 1; (25) 

(1 + lyo) + x{l - i/q) 

de-"' +x{d+l-de-'') 

51 x) = — , d>0. (26) 

d + x 

We fix i^o = 1/2 and I'l = d = 2. Here, d is chosen such that the introduced Minkowskian pole is located in the IR 
region. Thus, we have a function 

4"^('-)(Q^).f (27) 

which depends on the /3-function coefficients through fc„m and on A^. In Figs. [T] and O ^^^-"^'"^''^ and the perturbative 

coupling are plotted for = 2, 3, and 4, in the MS-scheme and for A = 0.4 GeV. In Fig. [T] one can see the very 

low dependence of on A^, and hence on f3i {i > 2), in the IR region. Fig. [5]shows that in the perturbative 

region, for a fixed A^, the perturbative and the analytized couplings merge. Besides, the A^ and Pi dependence of 
j^(N){iogs) moderate in the intermediate region. 
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FIG. 3: The parametrized part of A\ , given in Eq. I|28|l . is plotted for r; = hi = 1, A = 0.4 GeV, and various values of h2- 



We adopt the following point of view: the first term in Eq. (|23p gives a contribution to Ai (Q^) in the IR region 
with no free parameters (other than A), the only free parameters of the coupling being those contained in the second 
term of Eq. (|23l) .^ 

Now we turn to the second term of Eq. (|23p . Its parametric form is chosen without a physical motivation in mind. 
One just aims at a simple way of describing the coupling behavior in the IR. A possible choice is: 



e-^/(-) = '^i(^^|^e-^> (28) 

with three parameters. The parameter rj regulates its exponential suppression, hi is a global factor and gives the 
value of this term at a: = 0, and /12 describes, to a large degree, its behavior at x ^ 1. The effect of varying 77 can be 
observed by comparing Figures [3] and [H In each of these figures, the curves for various values of are plotted. 

In Fig. [51 curves of Af"^ for three different choices of parameters for are shown, illustrating various possible IR 
behaviors of the analytic coupling. 

The RS dependence of A^^"^ is contained only in the coefhcients knm, the parametrized contribution, Eq. (|28p . 
being RS independent. 



C. Minkowskian Coupling 



In this Subsection we present the Minkowskian coupling 2li(s) corresponding to the exponentially modified one- loop 
coupling. The Minkowskian (time-like) coupling 2li(s), for s > is defined from the Euclidean (space- like) coupling 
AiiQ^) by 



^ Of course, if we had made a different choice for Lo{x) and/or Li{x), the values of the free parameters in e f{x) would change. 
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FIG. 4: Same as in Fig. |3] with r] = 0.3 and hi = 1. 



-1 r"^'" dz 

2liW = / -Ar{-z), (29) 
where the integration path avoids the time-like semiaxis z > 0. The inverse relation is given by 



(30) 



The couplings 2li(s) and Ai(Cp) are related to each other in the same manner as the e+e" ratio R{s) and the 
Adler function D{Q'^) are (25l.l2^. The coupling Ai{Q^) is analytic in the complex plane excluding the Minkowskian 
semiaxis and behaves appropriately at infinity in order to have the spectral representation 



TT 7o cr + Q 



(31) 



where the spectral function pi = lm[Ai{—<J — ieA^)], with e ^ 0. Replacing Eq. ([5T|) in Eq. (I^^)) the following 
representation for the Minkowskian coupling is obtained 



2li(.) = i r-pi(a) 



(32) 



Returning to the evaluation of observables, the one-chain resummation integral in Eq. ([T]), which involves Ai in the 
integrand, can be written in terms of 2ti j2j] 



Jo 



(33) 
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FIG. 5: Three examples o{ Al ' , Eq. l|23ll, are plotted together with the perturbative coupling Opj' . The function a^^' is plotted only for 
q2 y ji^2 because it has a branch cut along < A^. The MS-scheme and A = 0.4 GeV are used. 



where F^{t) is the observable-dependent Minkowskian characteristic function. 

In this Subsection we present the spectral function and the Minkowskian coupling corresponding to the exponentially 
modified one-loop coupling 



1 / 1 



2x 



hi 



1 + X 



2-^ p-vV^ 



(34) 



Aogx 1-x {1 + i^o) + x{l - i^o) ^ (l + a;/2)2 

with X — Q^/A^ and vq — 1/2. The free parameters of the model are r], hi, and /12. The spectral function is 
divided in two terms, the A'^ = 1 contribution and the parametrized part of the coupling (proportional to hi): 
pi{a)^p[''-'\a)+p[^^'^^\a),with 



and 



Pi (o-) 



/?o(loga/A2)2 + 7r2 
2ae'''> 



/?o(a + A2) 

-I- sin(i'oA/o-/A)Re(^— 



(i^o%/^/A) 7r5((l + lyo) - (1 - t^o)CT/A2) 
1 

(l-hi/o)-(l-i^o)o-/A2 + ie 



(35) 



^(para)^^^ = 4/li(l - /isCj/A^) L ,5' ((j/A^ - 2) cos(77 V^/ A) -|- Rc f — 

L V (a 



(cr/A2-2 + ie)^ 



sin(77^/CT/A) 



(36) 



The Minkowskian coupling 2li(s) = 21^^ ^\s) + ^'f^'''^\s) is obtained by replacing the spectral function in Eq. 
The A'^ = 1 contribution is given by 
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1 1 



for s < (i±^)A2, and 



/5o7r Jo 



dx X sin(j/ox) 



2 _ i+£o a;2 + 1 



(37) 



aj^-^' (s) = ^ I + _ arctan 

PO ^ 2 TT TT 



1 1 



log(s/A2)N 2e 



a; sai{vQx) 



(38) 



for s > (y3^) A . Both expressions are logarithmically divergent as s ^ ( ) ^ ■ The parametrized contribution 
is given by 



2l(P"'^^)(s) = /ii 



2/ii 



dx X sin(?7a;) 



r 1 



4(7 



(a;2 - 2)2 



(39) 



for s < 2A2, and 



2((para)^^^ 



2/ii 



dx X sin(77x) 



(40) 



for s > 2A2. aJP'"''\s) behaves like (s - 2A2)-i for s w 2A2. 

In Fig. [H] the Minkowskian coupling 2ti is plotted together with the Euclidean coupling Ai for a fixed set of 
parameters. In the IR the both couplings differ drastically, while for energies ~ 2 — 3 GeV and higher the difference 
is just the usual (perturbative) 7r2-term [25l[26|. 

The Minkowskian coupling can be obtained also for higher loop values N . The resulting expressions are rather long 
and hence not shown here. For resummations, it is usually easier to use the Euclidean coupling. 



IV. EVALUATION OF OBSERVABLES 

A. Loop-level Dependence and RS Fixing 

For an energy scale relatively near to the boundary of the applicability region of pQCD, let us say Q — 4: GeV, 
'^pt~^^ and Op^"'*'' differ by ~ 2%, while the difference between and A[^^ is much smaller (see Figures [H and [S]) . 
Thus, in the perturbative region A^^"^ mimics Op^^ and the well established methods of pQCD can be applied using 
the analytic coupling A[^^ . 

As the value of Q decreases, the standard perturbative predictions for an observable ViCp) become unstable against 
variations in the number of loops and against the RS choice, and show a bad apparent convergence of the first terms of 
the series. All this is related to the strong dependence of the low-scale perturbative coupling under scale and scheme 
variations, and consequently also under the variation of N . In turn, this is related to the vicinity of the Landau 
singularities. For Q = 1 GeV the A^-dependence of the coupling can be seen in Fig. [1] for the MS-scheme. The 
couplings flp^"'^^ and Op^""'^ differ by « 40%. In contrast, the difference between the analytical couplings A^i^~^'^ 

and A^i^~'^^ is ~ 5%. In general, we shall call a coupling A^-{^^ "in the intermediate region N -independent coupling" 
if it varies at Q « 1 GeV by 5% or less when varying N (with > 2). 

Let us evaluate a space-like observable T) using truncated perturbation series at three- and four-loop-level, and 
using the analytic coupling A^i^\ In this Subsection we analytize the non- leading terms of the series making the 

replacement: (flp^')'^ [A^-^'^Y] however, a different approach will be discussed in the next Subsection. The three- 
and four-loop predictions for the leading twist contributions are: 



p(N=3) ^ aT +dMf'^? + d2[AT)\ 

P(N=4) ^ AfUdM'f^?+d2{Af^f+dMf^)\ 



(41) 

(42) 
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FIG. 6: The Minkowskian coupling Ul ' (s) = 2tf '(s) and Euclidean coupling Al '{Q'^) are plotted as a function of the energy scales y/s 
and Q, respectively. The case = 1 is considered, the relevant expressions are given in Eq. I|34| l and Eqs. I|37|l - H40| l. The parameters are 
fixed as follows; r] = 0.3, h\ = 0.1, and h2 = 0. The value A = 0.4 GeV is used. 



where, of course, I?^'^) and Ai are functions of Q^. In the intermediate energy region (Q ^ 1 GeV), where Ai 
can be of order one, powers of Ai (in particular the term d3(^^^-')*) are not necessarily suppressed with respect to 
A'f \ Hence, given the experimental value of V{Q^), the extracted values of A^^^ will depend strongly on N . This is 
problematic, because it is very unlikely to have TV-stable observable predictions together with a iV-unstable coupling. 
A possible way to obtain an in the intermediate re gion A^-independent coupling is to choose the RS separately for each 
observable (as in the effective charge method (27l l28l [2^). demanding for each observable: da = 0. Thus, extracting 

the coupling from a measured V, we obtain from Eq. (021) -^i'^^ = -^i''^ • On the other hand, from Eq. (^5)) 

_^(^.4) _ ^(A.=3) ^ ( ^ fc^^ log™ , J_. (43) 

m=0 " 

The latter difference is normally smaller than higher order uncertainties and experimental errors. We conclude that 
an in the intermediate region A^-independent coupling A^-^\ used together with an observable-dependent RS choice, 
provides a A^-stable framework for the evaluation of observables at intermediate scales which is consistent within the 
errors. 

Let us shortly review how within the evaluation method of Milton, Solovtsov, Solovtsova and Shirkov (MSSSh) 
P, 0, 0, 01 the A^-stability problem of observables at intermediates energies is solved. In this approach, an observable 
in a given RS at A^-loop-level is expressed as: 

WN) _ .(AT) , .(A) , AN) 

'^MA - A, MA + "1-^2, MA + ■ ■ • + "W-l-^jv,MA- \^^) 

The "minimal analytic" (MA) functions A'^^^^iQ'^) [E [3, H, H are computed using dispersion relations, where the 
spectral functions are obtained from the functions {a'^\Q'^))^ . Hence, perturbative and MA couplings possess the 
same discontinuity along the Minkowskian semiaxis. The behavior of A^j^-^pJ^Q'^) for large is the same as the one 
of 

('ipt^''(*9^))'' j modulo (unwished) power-behaved terms. The analytization procedure is unique, in the sense that it 
does not introduce any new parameter. Thus, the IR behavior of the couplings is not modeled but given by the first 
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coefficients of the /3-function; from our point of view, this a limitation of the approach. An interesting feature of the 
MSSSh procedure, connected with the TV-dependence of the couphng, is that the couplings -4^,^^ (with k > 2) are 

suppressed compared to (^i^a)'^- ^ consequence, the analytized truncated perturbation series is rather stable 
against variations in the number of loops and against the RS choice. It also shows well the apparent convergence of 
the first terms of the series. 



B. Skeleton-motivated Approach 



In this Subsection we present the skeleton- motivated approach for the evaluation of observables [llj, [l^] . In order 
to be specific, we consider the next-to-next-to-leading order case, other approximations being generalizations of it. 
The third order renormalization-group (RG) improved truncated perturbation series has the form: 

2?pt(Q') = apt(Q') + di a^t(Q') + ^2 a^t(Q'), (45) 
where the perturbative coupling Opt obeys the third-order RG equation: 



da 



pt 



91ogQ2 



(46) 



The skeleton expansion exists in QED if light-by-light subdiagrams are excluded [30|. In QCD, its existence is not 
certain, but it can be postulated [3ll.l33|. The evaluation approach is motivated by the skeleton expansion but it does 
not rely on its existence. The skeleton expansion is given by: 



/>00 7 1 OQ /"OO Ji 



(47) 



where the Euclidean characteristic functions F^{ti,. . ., i„) are symmetric functions and are normalized as: 



f'^FUt)^!, J'^^FUtuh)^!, (48) 

and sf are the skeleton coefficients. It is assumed that the characteristic functions do not have any dependence on 
the running coupling constant. If the representation (|47p exists, its expansion in powers of apt{Q'^) must correspond 
to the all order perturbation series of the considered observable. The IR contributions of the integrals in Eq. ([T7|) 
suffer, as we discussed earlier, from ambiguities due to the singularity and cuts of the perturbative coupling constant. 

The analytization is formally made by replacing the perturbative coupling apt{Q'^) of all skeleton integrands in 
Eq. (l47l) . by an analytical coupling Ai{Q^). The replacement fixes the IR ambiguities of the integrals. In order to 
obtain the analytic truncated series analogous to Eq. (|l5)) . the function Ai{Q^) is Taylor-expanded around (here, 
for simplicity) t = 1 inside all skeleton terms. In the considered third order treatment we need the first and second 
derivatives of Ai{Q^). Instead of the derivatives we use equivalently the analytical functions A2{Q^) and A3{Q^), 
which are defined by: 



dAi 
91ogQ2 

d^Ai 

5(l0gQ2)2 



= -[/3o^2(Q')+/3i^3(Q')], (49) 
= 2/32 ^3 (Q'). (50) 



The latter equations are the truncated RG equation and the first derivative of it (in higher order treatments, higher 
derivatives are involved), with the replacement ap^. 1-^ Ak- Perturbatively, the differences A2 — Opt and .43 — ap^. are 
of O(apj). In Fig. [71 the couplings A2 and A3 obtained using the analytic coupling without power corrections defined 
in Subsection IIII Bl are plotted . 

After Taylor-expanding, the terms corresponding to order a^^ (e.g. A2 or ^1^3) and higher are dropped and the 
skeleton integrals are formally performed. As a result, we obtain: 
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FIG. 7: The couplings ^2 and .43, defined in Eqs. II49II and Il50|l . together with the corresponding couphng Ai, defined in Subsection 
nil Bl are plotted as a function of Q, in the MS-scheme, with A = 0.4 GeV. The parameters used for the couplings are r} = 0.3, hi = 0.1, 
and /i2 = 0. 



2?an(Q') = Ai{Q^)+[f3ofr{l)A2{Q') + sfAj{Q^)] (51) 
+ [iPofFi^) + /3i/f (1))^3(Q') + 2sf(3ofr,oa)AiiQ^)A2{Q') + sfAliQ^)] , 



with the momenta 



f?{t*)= I jFUt)i-logt/U)\ 



/dtl dt2 „e , , , 



(52) 



Thus, comparing Eqs. (|45p and (|5ip. we see that in the analytization process apt is replaced by Ai, ap^. is replaced 
by a linear combination of A2 and Ai, and a^^ is replaced by a linear combination of ^3, A1A2, and Ai- In the UV 
regime, the difference between both equations is of order Op^. The coefficients and sf and the values of /f , 
and /^Q can be obtained in a given RS, for the third order treatment, from the perturbative coefficients di and ^2 
using their dependence on the number of flavors Nf (see [HI, [T2l|). 

The definitions (j49|) and ([50]) are made in such a way that, the result of the third order Taylor-expansion of the 
analytic coupling inside the one-chain term 

°° df 

- FUt)Ai{tQ') ^ Ai{Q') + f3ofTWA2{Q') + {Plf^{l) + Mr{l))A3{Q^), (53) 







corresponds to the perturbative third order truncated expression 

^ FUt) apt(ig^) = apt(Q2) + PofFaKiQ') + iP'of^i^) + /3i/f (l)X(Q') + 0{a%), (54) 
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with Opt I— > Ak- A different cfioice of A2 and A3 would not allow this correspondence. 

If the exponentially modified coupling of Section IIIII is used in Eq. (|5ip , no power-behaved terms are introduced. 

On the other hand, if the characteristic function Fj,{t) is known, the object we are interested in is the analytic 
version of the one-chain resummed extension of Eq. (|45p . In order to obtain it we expand the skeleton integrals as 
before, except for the leading term which we keep unexpanded. We obtain 

r°° rlf 

V„{Q^) = / - FUt) AiitQ^) + sfAliQ^) + [s^AliQ^) + 2sf /3ofToi^)Ai{Q')A2{Q')] . (55) 

Besides, as we motivated in the previous Subsection, in order to decrease the TV-dependence of the coupling, we adopt 
an observable-dependent RS choice, such that each observable is written as 

T^iQ') = / - FUt) AiitQ^) + sfAliQ^) + On, (56) 
Jo ^ 

i.e., evaluation in a RS where and f^Q are equal to zero (again, see [Ulil for details). The error, for the present 
(third order) case, is formally of order four. On = O4. The final expression Eq. (|56p remains unchanged in higher 
order treatments, the error 0„ being one order of magnitude higher than the order of t he p erturbative expression. 

This is one variant of the evaluation approach. Other possibilities can be found in [lll.lT3|. In particular, a different 
RS can be chosen [l^] such that the second term of Eq. ([5S|) is replaced by the first derivative of the analytic coupling 

-D'iQ') = fj Fiit) Mm + + On, (57) 

where the derivative of Ai can be written in terms of A2 and ^3 according to Eqs. and ([501) . Refs. [ill, 
this evaluation approach is applied to the Adler function, the Bjorken polarized sum rule and the semihadronic r 
decay ratio, using three different models for Ai (which do introduce running coupling power-terms). 



V. CONCLUSION 



Motivated by the belief that relevant non-perturbative QCD contributions arise only from the integration of 
IR degrees of freedom, a new class of models for an analytic QCD coupling is proposed. The main characteristic 
of the proposed analytic coupling is that its difference to the perturbative one, in the UV region, is smaller than 
any inverse power of the energy. As a consequence, a general analytization procedure of observables which makes 
use of this coupling yields predictions which have no power-behaved contributions arising from UV degrees of 
freedom. In addition, the constructed coupling is finite when Q — > 0. The analytic coupling is conceived as the 
central object in the construction of observables with the required analytic properties. A particular procedure is 
presented in Subsection lIVB| the skeleton- motivated approach for the evaluation of observables, developed in [lllIT^. 

For a coupling whose difference to the perturbative coupling is a very suppressed power term ^ (A^/Q^)''™"" 
(fcmax ^ 1), it is shown that in the one-chain resummation of space- like observables 'D{Q^), the UV regime still 
contributes an appreciable power correction ~ (A^/Q^)". Here z = n is the location of the leading IR renormalon 
of V in the Borel plane, and therefore, the leading OPE contribution to this observable is also ~ (A^/Q^)". Thus, 
even the suppressed difference 5Ai{Q^) ~ (A^/Q^)*^"""' cannot be accommodated within the ITEP interpretation 
of the OPE, which says that the power-suppressed terms in T>{Q^) are contributions from the IR regime only. We 
emphasize that the last conclusion is obtained using the criterium for the identification of UV contributions, in the 
one-chain resummation integral, described in Section [ill It is based on comparing the momentum flowing through 
the photon/gluon internal propagator, with the characteristic scale of the strong interactions, A. 

An analytization procedure that makes use of the proposed coupling, must be considered as a phenomenological 
attempt with the main motivation to enlarge the energy range of applicability, and improve the description capability 
and the predictive power of perturbative QCD. Of course, we are most interested in the intermediate energy region 
Q fa 1 — 2 GeV. The model includes, in addition to the QCD scale, a number of parameters (three parameters) 
describing the IR behavior of the coupling. In the intermediate energy region, the proposed coupling has low N 
(loop-level) and RS dependence. By comparing predictions with experimental data for low-energy observables, 
the three parameters can in principle be fixed and the low-energy behavior of the coupling can be obtained. If, 
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in addition, an all order resummation is made, as in Eqs. (|56p or (|57|) . then values of the coupling at low Q are "probed". 

In Subsection IIV Al we motivate an observable-dependent RS-fixing criterium. It provides a consistent A^-stable 
framework for the evaluation of observables, using a coupling with low A^-dependcnce at intermediate scales. 

The evaluation approach given by Eq. ([56]) or by Eq. ((57|l generates power-behaved terms only of IR nature. 
However, it is likely that the model cannot describe all non-perturbative effects relevant to low-energy observables. 
For example, chiral symmetry breaking effects are probably not accounted for within this framework. This motivates a 
mixed approach. Namely, the one where to Eq. ([56)l or ((57| OPE power-behaved terms are added. As in the standard 
one-chain resummation, the expectation is that non-perturbative contributions generated by the analytic coupling 
correspond to an important part of the non-perturbative effects otherwise carried by OPE terms in the standard 
(perturbative QCD plus OPE) approach. 
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